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ABSTRACT 

We present results of the multicolour UBVR photometry of the high-amplitude EC 14026- 
type star, Balloon 090100001. The data span over a month and consist of more than a hun- 
dred hours of observations. Fourier analysis of these data led us to the detection of at least 
30 modes of pulsation of which 22 are independent. The frequencies of 13 detected modes 
group in three narrow ranges, around 2.8, 3.8 and 4.7 mHz, where the radial fundamental 
mode, the first and second overtones are likely to occur. Surprisingly, we also detect 9 in- 
dependent modes in the low-frequency domain, between 0.15 and 0.4 mHz. These modes 
are typical for pulsations found in PG 1716+426-type stars, discovered recently among cool 
B-type subdwarfs. The modes found in these stars are attributed to the high-order g modes. 
As both kinds of pulsations are observed in Balloon 090 100001, it represents a link between 
the two classes of pulsating hot subdwarfs. At present, it is probably the most suitable target 
for testing evolutionary scenarios and internal constitution models of these stars by means of 
asteroseismology. 

Three of the modes we discovered form an equidistant frequency triplet which can be 
explained by invoking rotational splitting of an £ = 1 mode. The splitting amounts to about 
1.58 fiHz, leading to a rotation period of 7.1 ±0.1 days. 
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1 INTRODUCTION 

Soon after the discovery of about a dozen pulsating B-type sub- 
dwarfs (sdB) by the South Afric an Astronomical Observatory as- 
tronomers I Kilke nny et ail i 1997b and the next papers of the series], 
two other groups (Billeres et al. 2002, 0stensen et al. 2001a,b, 
Silvotti et al. 2002a) undertook extensive searches for these in- 
teresting multiperiodic pulsators, called — after the prototype — the 
EC 14026 stars. The searches resulted in finding further variables, 
so that there are 32 EC 14026 stars known at present. Their pul- 
sation periods range between 1.5 and 10 minutes, with a typical 
value of 2-3 minutes. The amplitudes are rather small and rarely 
exceed 10 mmag. Some extreme cases are known, however. The 
longest periods are observed among the members having the low- 
est surface gravity: V338 Ser = PG 1605+072 (Koen et al. 1998a, 
Kilkenny et al. 1999), KL UMa = Feige 48 (Koen et al. 1998b, Reed 
et al. 2004a), HK Cnc = PG 0856+121 (Piccioni et al. 2000, Ulla et 
al. 2001), V2214 Cyg = KPD 1930+2752 (Billeres et al. 2000), 
V391 Peg = HS 2201+2610 (Ostensen et al. 2001a, Silvotti et 
al. 2002b), and HS 0702+6043 (Dreizler et al. 2002). Some of them 
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(e.g., V 338 Ser arid HS 702+6043) have large amplitudes. Re- 
cently, Oreiro et al. (2004) discovered another high-amplitude and 
coolEC14026 star, Balloon 090100001 = GSC 02248-01751 (here- 
after Bal09). The semi-amplitude of one of the detected modes 
amounts to 60 mmag in white light, the largest amplitude known 
in an EC 14026 star. With B ~ 1 1.8 mag, the star is actually one of 
the brightest members of the group. 

B-type subdwarfs are believed to be low-mass (~0.5 Mq) 
evolved stars that burn helium in their cores (Heber 1986, Saffer 
et al. 1994). In the H-R diagram they populate the extended hori- 
zontal branch (EHB) region. As their hydrogen envelopes are very 
thin, they do not ascend the asymptotic giant branch and do not pass 
the planetary nebula phase of evolution iGreggio & Renzini 1990). 
Instead, they evolve almost directly to the w hite-dwarf stage, pass- 
ing only through the hotter sdO phase iDorman et all 1 995). Their 
evolutionary past, however, is not well understood. The discovery 
of pulsations in sdB stars opened therefore a new way for investi- 
gating their internal structure by means of asteroseismology (see, 
e.g., Brassard et al. 2001). The effective temperatures of EC 14026 
stars range from 29 000 to 36 000 K, their surface gravities, log g 
= 5.2-6.1. They do not occupy, however, a well-defined instability 
strip in the log T c ff-log g plane but are spread among non-pulsating 
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sdB stars. Curiously enough, their pulsations were predicted almost 
at the time of discovery (Charpi n'et et alJ 19961) as corresponding to 
radial and nonradial low-order low-degree p modes driven by the 
k mechanism. The driving comes from the metal opacity bump 
that originates from the ionization of heavy metals, mainly iron 

(Charpinet et al. 1997, 2001); 

About two years ago, iGreen et"al] 120031) discovered low- 
amplitude multiperiodic light variations with periods of the order 
of one hour among the coolest (T e g < 30 000 K) sdB stars. The 
periods of the members of this new class of variable stars indi- 
cate that they are high-order g modes. Although lGreen et aljfcoojl) 
announced the discovery of 20 stars of this type, detailed anal- 
ysis was so far p ublished only for the prototype, PG 1716+426 
iReed et alj 2004a). The pulsations were explained by means of the 
same pulsation mechanism as the short-period ones l Fontai ne" et alJ 
2003), but with degrees t = 3 or 4. As far as the effective temper- 
ature and the character of excited modes are concerned, the two 
classes of variable sdB stars resemble younger main-sequence pul- 
sators, /3 Cephei and SPB stars (see, e.g., Pamyatnykh 1999). In the 
log T c ft-log # plane, the areas where the EC 14026 stars and their 
long-period counterparts are found slightly overlap. It seems there- 
fore possible that both kinds of pulsations observed in sdB stars 
might coexist in a single star of this type. 

The sdB nature of Bal09 was discovered by iBixler et alJ 

ll99ll) using the balloon-borne SCAP ultraviolet tel e scope . Its 
short-period variability was detected by Oreir o et"all 120041) by 
means of the 80-cm IAC80 telescope at Observatorio del Teide 
(Tenerife). The effective temperature of 29 500 K and log g ~ 5.3, 
derived by the same authors, place it among the coolest and most 
evolved EC 14 026 stars. Fro m the data obtained on a single obser- 
vational night. lOreiro e t al. 1 2004) were able to detect only two in- 
dependent modes. Owing to the large amplitude and relatively long 
periods, the star was an obvious target for follow-up time-series 
photometry. We therefore decided to carry out multicolour obser- 
vations of this star. New observations were also obtained by the dis- 
coverers (Oreiro, private communication). In this paper, we present 
the results of our photometric run, covering over a month. We show 
that Bal09 exhibits both kinds of pulsations discovered in sdB stars. 



2 OBSERVATIONS AND REDUCTIONS 

We started observing Bal09 on June 30, 2004. The bulk of the 
data was, however, obtained between August 17 and Septem- 
ber 19, 2004, on 18 nights. All these observations were obtained 
at Mt. Suhora Observatory with a 60-cm reflecting telescope, 
equipped with the ST-10XME CCD camera and Johnson-Cousins 
U BVR filters. The camera field of view was 6.8 x 4.5 arcmin and, 
in order to minimize storage time, the 3x3 pixels on-chip bin- 
ning was employed. Except for the first night, June 30, when only 
R filter was used, and the second night, when V-filter measure- 
ments were not made, the observations were carried out through 
four filters, U, B, V, and R. The mean exposure times amounted 
respectively to 15, 10, 7 and 7 s, resulting in about 1 minute-long 
filter cycle. This enabled us to obtain roughly 6 points per main 
pulsation period of Bal09 in all four filters. In total, 125 hours of 
observations were obtained in the U and B filters, 119 in V, and 
126 in R. 

Additional observations were carried out on three observing 
nights between September 22 and 25, 2004, at Loiano Observatory 
with the 1.52-m Cassini telescope and the BFOSC CCD camera 
covering 13 x 12.6 arcmin on the sky. In total, we gathered 13 
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Figure 1. Differential U BVR magnitudes of Bal09 obtained during our 
run, plotted against heliocentric Julian Day. The observations were obtained 
at Mt. Suhora Observatory, except those on the last three nights in V that 
were taken at Loiano Observatory. 



hours of observations only in filter V. The exposure times were 10 
or 15 s long, depending on sky conditions. The time distribution of 
all observations is shown in Fig.Q 

Over 7000 individual frames were obtained in U, B and R 
and about 8500 in V. All these frames were corrected for bias, 
dark and Hatfield in a standard way. Then, the magnitudes of all 
stars detected in a frame were obtained by means of profile-fitting 
procedures of the Daophot package 1 Stetsonl ll987l) . The aperture 
magnitudes were also derived. In order to minimize photometric 
errors, the apertures were scaled with seeing. As a measure of see- 
ing, we used the <rs of the two-dimensional Gaussian function that 
described the analytical part of the point-spread function (PSF) cal- 
culated by Daophot. The positions of the stars were also taken from 
the PSF fits. The seeing scaling factor, optimal from the point of 
view of the photometric errors, depended on the aperture magni- 
tude of a star. We derived the character of this dependence from 
the post-fit residuals on a single good night and then applied it to 
all observations. As a comparison star, we chose the nearby, rel- 
atively bright star, GSC 02248-00063. It is brighter than Bal09 in 
B, V, and R, but slightly fainter in U. As can be seen in Fig.Q 
the mean differential magnitudes range from —0.23 mag in U to 
1.54 mag in R. This large range reflects large differences in colour 
between Bal09 and the comparison star. Surprisingly, there is prac- 
tically no difference between the mean differential magnitudes for 
the V-filter data from Mt. Suhora and Loiano. This means that the 
instrumental systems in both sites are very similar and the V data 
can be safely combined. 

The differential aperture photometry showed less scatter than 
the PSF photometry. Consequently, for the frequency analysis that 
follows, we used the aperture photometry. One night's sample dif- 
ferential magnitudes in V are shown in Fig.|2| 



3 FREQUENCY ANALYSIS 

Prior to the frequency analysis, we rejected data points that had 
photometric errors larger than a certain, arbitrarily chosen value, 
viz., 0.04 mag for U, 0.03 mag for B, and 0.025 mag for V and R. 
As the amplitudes of the light variation of Bal09 are quite large, the 
next three steps, i.e., rejecting outliers, correcting for second-order 
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Figure 2. Observations of Bal09 on a single night (September 6, 2004) ob- 
tained at Mt. Suhora Observatory in filter V . The solid line is the synthetic 
light-curve consisting of 30 periodic terms listed in Table 1 . 



extinction effects and detrending, were done using the residuals 
from a preliminary fit. For this purpose, sequential prewhitening 
was done, until all periodic terms with signal-to-noise (S/N) ratio 1 
larger than 6 were subtracted from the data. 

Having obtained the residuals, we first rejected all outliers 
from the original data employing 5-ct clipping. Since the compar- 
ison star is much redder than Bal09, and we made observations 
through wide-band filters, second-order extinction effects can be 
expected. In order to correct for them, we plotted the residuals, 
Am les , obtained earlier, against the air mass, X, and derived the 
coefficients of the relation 



Am r , 



a(X - 1.25) 



by means of the least-squares method. The linear coefficient, a, was 
significantly different from zero and amounted to +0.044, +0.033, 
+0.009, and +0.012 ± 0.001 for U, B, V, and R, respectively. The 
a(X — 1.25) factor was then subtracted from the original data. Fi- 



We define S/N as the ratio of the amplitude of a given peak in the Fourier 
periodogram to the average amplitude in the whole range of the calculated 
frequencies. This definition differs from the common definition of this value 
in the sense that usually the noise N is calculated from the periodogram 
of the very last residuals. As we use S/N only for deciding when to stop 
prewhitening, this difference in defining S/N has no influence on the final 
result. 
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Figure 3. Fourier amplitude periodogram of the detrended V-filter data of 
Bal09. The highest peak at about 2.8 mHz is truncated; its height is equal 
to about 53 mmag. 



nally, the prewhitening was repeated and new residuals were com- 
puted. 

Because 90% of our data in V and all data in the remaining 
three filters come from the same telescope/detector combination, 
they constitute a very homogeneous dataset. Nevertheless, some in- 
strumental effects can produce artificial trends in the data, increas- 
ing amplitudes at the lowest frequencies, making it difficult to say 
which low-frequency signals are real. We will discuss this prob- 
lem in Sect. 3.7 using data that were not detrended. At the present 
stage, detrending was performed, again using the residuals. First, 
the average magnitudes in 0.07-day intervals were calculated and 
then interpolated using the cubic spline fit. This fit was subtracted 
from the original data. The 0.07-day interval was chosen by trial 
and error. This choice secures that all signals at frequencies below 
~0. 1 mHz are effectively reduced while those with higher frequen- 
cies, where g modes are observed (Sect. 3.2), remain unaffected. 



3.1 Detected frequencies 

Detrended data were analyzed by means of the standard proce- 
dure that included calculating the Fourier amplitude periodograms 
and sequential prewhitening of detected signals. At each step of 
prewhitening, the parameters of sinusoidal terms found earlier were 
improved by means of the non-linear least-squares method. Data in 
each filter were analyzed independently. The periodogram of the V- 
filter data, up to the frequency of 10 mHz, is shown in Fig. [3] The 
structure of the aliases can be seen in Fig.|4| where we show Fourier 
periodogram of artificial data containing a single sinusoidal term 
with frequency and amplitude of the main mode and distributed 
in time as the real observations in V . The daily aliases are quite 
high, as both sites where the V observations were carried out are 
located at a similar longitude. Therefore, the frequencies we derive, 
especially the low-amplitude ones, may suffer from the 1 (sidereal 
day) -1 = 11.6 /iHz ambiguity. Our data cov er 35 days, re sultin 
in a f requency resolution of about 0.5 /iHz lLoumos & Deemin 
1978). 

The Fourier spectrum of Bal09 shown in Fig.[3]is dominated 
by the strong peak around / ~ 2.8 mHz and its daily aliases. It 
is the main frequency discovered bv lOreiro et alj i2004l) . Signals 
at the harmonics, 2/ (5.6 mHz) and 3/ (8.4 mHz), can also be 
seen. The second frequency detected by these authors, at about 
3.8 mHz, is also well seen. In addition, we see low-amplitude peaks 
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Figure 4. Fourier amplitude periodogram of the artificial data containing a 
single sinusoidal term with frequency of 2.8074594 mHz and amplitude of 
53.3 mmag, distributed in time as the real observations in V. 

at 4.7 mHz and, surprisingly, at low frequencies below 0.4 mHz. 
N ote that the peak at frequency of 0.3 mHz was already detected 
bv lOreiro et alj i2004t) . but these authors were not able to attribute 
it convincingly to the star because of the very short time interval of 
their observations. 

The analysis was carried out until all terms with S/N > 4.5 
were extracted. The number of periodic terms found in this way 
was not the same in all filters, mainly due to the differences in the 
detection threshold. The threshold was the lowest in V, where we 
found 27 independent terms and 7 combinations terms, including 
harmonics. In U, we found 20 + 3, in B, 22 + 3, and in R, 23 + 6 
terms. Of all these terms, 17 (including 3 combination terms) were 
found in all four filters. We finally adopted a model consisting of 22 
independent terms and 8 combination terms. The model comprised 
all terms detected in V and at least one other filter, and all com- 
bination terms. The frequencies of all terms included in this model 
and the results of the fit of this model to our multicolour time-series 
data are presented in TableQ 

It is clearly seen from Table 1 (see also Fig.0 that the frequen- 
cies of modes detected in Bal09 cluster in four narrow frequency 
ranges: below 0.4 mHz, around 2.8, 3.8, and 4.7 mHz. After remov- 
ing the 30 terms from the data, the Fourier spectra of the residuals 
still exhibit an excess power in these four ranges (Fig.|5J- Because 
the noise level increases towards low frequencies, we stopped ex- 
tracting terms with frequencies below 1 mHz. However, extracting 
all terms down to S/N = 4.5 for frequencies in the range 1-10 mHz, 
led us to detecting additional 9 low-amplitude terms. Their frequen- 
cies are listed in Table [2] Note that only some are found in two 
bands, and even in these cases there is an ambiguity in frequency 
because of aliasing. It is clear that multi-site data are necessary to 
verify the reality of these terms. 

We will now discuss frequency contents of the multicolour 
data of Bal09 in the four frequency ranges mentioned above, illus- 
trating the discussion with the Fourier amplitude spectra in the V 
filter obtained at different steps of prewhitening. 

3.2 Low-frequency region (0.15-0.4 mHz) 

Out of 22 independent terms included in the model listed in Ta- 
ble 1, nine were found in the low-frequency range, between 0.15 
and 0.4 mHz. Five (/a to /d and /i) were detected in all four fil- 
ters. The Fourier spectrum in this range (Fig.|6j shows complicated 
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Figure 5. Fourier spectra of the UBVR residuals after removing 30 terms 
listed in Table 1 . 



Table 2. Additional high-frequency terms with S/N > 4.5 found in the resid- 
uals. 
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structure; the first four terms have amplitudes of about 2 mmag 
in V, the remaining ones are even weaker. The dense spectrum of 
terms with low amplitudes in this range causes that aliasing is really 
a problem here. 

3.3 The 2.8 mHz region 

The strong peak at frequency of about 2.81 mHz, found by 
lOreiro et alj 120041) . is resolved in our data into several compo- 
nents. The strongest component has frequency /i = 2.8074594 ± 
0.0000006 mHz (period 356.2 s) and is non-sinusoidal in shape, 
as we detected two its harmonics. All combination frequencies we 
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Table 1. Frequencies, UBVR semi-amplitudes, and phases of 30 sinusoidal terms adopted in the pulsation model of Bal09. Terms in a given frequency 
range are listed according to decreasing V amplitude. Amplitudes in square brackets indicate that the term was not detected in the data obtained in a given 
filter. Phases are given for the epoch HJD 2453252.0, the errors of their last digits are given in parentheses. The entries at the bottom of the table give global 
parameters of the fit. N a b s denotes the number of the data points in a given filter, cta , the rms error of the semi-amplitude, and Am, the mean differential 
magnitude (variable — comparison). Detection threshold corresponds to S/N = 4.5 in the Fourier periodogram of the residuals. 
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detected, except + fs, involve /i. The 6-minute oscillations cor- 
responding to fi dominate the light curve ofBal09 (see Fig.|2j. We 
detected six terms in the vicinity of /i that were included in the 
30-term model. One of the most interesting features seen in Fig.0 
is that /2, /3, and form an equidistant frequency triplet. The 
frequency differences between the components of the triplet are the 
same within the errors and amount to: fs — / 2 = 1 .574 ± 0.003 /iHz 
and / 4 — / 3 = 1.585 ± 0.008 /jHz. Thus, they can be resolved only 
if observations cover more than a week. It is natural to suggest that 
the triplet is a rotationally split one; we will return to this inter- 
pretation in Sect. 5.3. The beat period between /i and the triplet 



components is of the order of 0.7 d, and this beating is responsible 
for the changes of amplitude in the observed light curve (Fig.[2j- 

At a frequency about 0.03 mHz larger than that of the triplet, 
three next frequencies were resolved. They are closely, although 
not equally spaced, f 7 - f 5 = 1.72 ± 0.04 ^Hz, f 6 -f 7 = 2.85 ± 
0.04 /^Hz. 

One additional term was found in the 2.8 mHz region in the 
residuals of the 30-term model (Table|2j- This /14 term, found only 
in the V data (Fig.0, lies very close to the combination frequency 
fi + f-i — fs sa fi + fz — /4, but if real, it is more likely to be an 
independent mode. The reason is that a much higher amplitude is 
found when solving for this frequency rather than assuming that it 



6 A. Baran, A. Pigulski, D. Koziel, W. Ogloza, R. Silvotti and S. Zola 




Figure 6. Fourier amplitude spectra of the V data of Bal09 in the low- 
frequency domain. The number on the right shows how many periodic terms 
were extracted from the data before calculating the periodogram. Vertical 
dashed lines correspond to frequencies that have been already extracted. 
The inverted filled triangles indicate frequencies included in the model 
given in Table 1 , the open triangle, that detected only in V data. 



has a combination value. The difference /i 
± 0.03 fiHz. 



/i4 is equal to 1 .46 



3.4 The 3.8 mHz region 

A peak at 3.78 mHz was already seen by lOreiro et alj |2004) in 
their Bal09 observations. We resolve it into three terms, counting 
only those that appear in TableQor nine, if those from Tableware 
taken into account. It is interesting to note that /20 — /19 = 1.61 
± 0.04 /iHz, a value that, to within the errors, is the same as the 
separation of the three frequencies around 2.825 mHz. As can be 
seen from Fig. [H] (but see also Fig. [5}, the region is rich in low- 
amplitude modes. Not all can be identified unambiguously from 
our data due to aliasing and the fact that many appear only in one 
periodogram. 



3.5 The 4.7 mHz region 

Mode s with freq uencies around 4.7 mHz were not detected by 
lOreiro et alj 120041) . We included three terms from this region in 
our 30-term model, but two more are detected with S/N > 4.5: one 
in B and one in V (Fig.|5J- All terms in this region have very small 
amplitudes, around 1 mmag or even smaller. 



3.6 Combination frequencies 

The combination frequencies appear in three distinct regions (see 
Fig. 1101 and Table0. The richest is the 5.6 mHz region where the 
2/i harmonic and the sums of the highest-amplitude terms occur. 
Five combination frequencies are detected in this region, including 
the sum of /1 and all three triplet components, as well as /2 + /s. 

In the 8.4 mHz region, the 3/i term is detected only in V. 
However, although not detected above S/N = 4.5, the 2/i + /2 is 
also very well seen in the periodogram of the residuals ("Fig. llQI . 

For the assessment of reality of the signals, the most impor- 
tant is, however, the detection of the difference term /1 — /b at 



30 
30 
10 



15 

I" 
5 

P 

5 



o 

E 



3 
2 
1 




in 
in 



in 
in 
in 



■9t 




-34 f 



III 



2.78 2.8 2 82 2.84 2.86 

Frequency [mHz] 



2.88 



Figure 7. Fourier amplitude spectra of the V data of Bal09 in the vicinity 
of the main mode, /1 = 2.80746 mHz. Notation is the same as in Fig.l6l 
except that the frequency listed in Tablel2lis indicated by the inverted gray 
triangle. Note the differences in the ordinate scale. 



2.441651 mHz (see Fig. llQI top panel). This combination term in- 
volves frequencies from the high and low-frequency ranges. Its oc- 
curence proves that both /1 and /b originate in the same star. In 
Fig.[l0| we have also indicated the difference combination terms 
between /1 and three other low-frequency terms with the largest 
amplitudes. A detailed look at the residual spectrum in V shows 
that in addition to the peak at /1 — /b there is also a low but still 
noticeable peak at /1 — fc = 2.48165 mHz. 



3.7 The lowest frequencies 

At the beginning of this section we mentioned that detrending was 
performed prior to the Fourier analysis. This caused all information 
on signals with frequencies below ~0.1 mHz to be lost. However, 
some real low-frequency signals may be present in the photome- 
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Figure 8. Fourier amplitude spectra of the V data of Bal09 around fre- 
quency 3.8 mHz, showing representative steps of prewhitening. Symbols 
are the same as in Figs.lo1andl7l Note that out of six additional frequen- 
cies (Table 0, indicated by the gray triangles, two were found in the V 
data only. These two are shown in the second panel from the bottom, the 
remaining four, in the bottom panel. Triangles indicating alias frequencies 
are connected with horizontal lines. 
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Figure 9. The same as Fig.|8] but around frequency 4.7 mHz. 
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Figure 10. The same as Fig.|8] but in three regions where the combina- 
tion frequencies occur. See text for the explanation of the difference terms 
indicated in the top panel. 
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Figure 11. The Fourier amplitude spectra of the Bal09 U BVR observa- 
tions for frequencies below 0. 1 mHz. The data were not detrended, but all 
terms listed in Tables[T1andl2lwere removed. 



try of Bal09. In particular, we can expect the difference combina- 
tion frequencies. Therefore, we decided to analyze also the data 
that were not detrended. First, the already known terms with fre- 
quencies higher than 0. 1 mHz were removed from the data, and the 
Fourier periodograms of the residuals were calculated. They are 
shown in Fig. ll ll for all four filters. 

We see that peaks close to and 1 d _1 = 11.6 ^iHz appear in 
all periodograms indicating that some night-to-night changes, pre- 
sumably of instrumental origin, are present in the data. This can be 
seen most clearly in the R data. However, there is an indication of 
a significant peak at frequency of about 9.26 /jHz or its alias, 11.57 
— 9.26 = 2.31 ^iHz in the U and B data. The corresponding peri- 
ods are equal to 1.25 and 5 days. It is difficult to say whether this 
is a real periodic change. We think it is rather an artefact caused by 
transparency changes and/or instrumental effects. The differential 
combination terms are not visible below 0. 1 mHz; they are proba- 
bly lost in the noise. 

3.8 Reality of the detected modes 

It is very important from the point of view of the results of this 
paper to be certain that the modes we detected, especially the low- 
frequency ones, are real. Moreover, we would like to be certain that 
they originate in Bal09, and not in the comparison star. In order to 
verify the latter, we analyzed differential photometry of the com- 
parison star with respect to another star in the field, GSC 02248- 
00425, called thereafter the check star. Unfortunately, this star is al- 
most 4 mag fainter than the comparison star, so that only the V and 
R photometry is usable for it. The periodograms of the differential 
photometry of the comparison star with respect to the check star 



show no peak above S/N = 3.7 corresponding to the semi-amplitude 
of 3.0 mmag in V data and no peaks above S/N = 4.1 correspond- 
ing to 1.9 mmag in the R data. This justifies that all modes with 
semi-amplitudes exceeding 2 mmag in R or 3 mmag in V originate 
indeed in Bal09. Among them, there are two low-frequency modes, 
/a and /b . We cannot, however, exclude the possibility that some 
low-frequency modes originate in the comparison star. 

An additional argument for the reality of the low-frequency 
modes may come from an independent detection of a given mode 
in the data split into two parts; see lKilkennv et alii 19991) for an ex- 
ample. Consequently, we divided the photometry of Bal09 into two, 
roughly equal parts, before and after HJD 2453254. Then, the anal- 
ysis of the two parts was performed independently. This was done 
for data in all four filters. Of 22 independent terms from TableQ we 
detected 14 terms in both parts with S/N > 4.5, or 16 if detection 
threshold was lowered to S/N = 4. The remaining six terms, /f, /g, 
fe, fi, fn, and /12, were detected in the periodogram of only one 
part of the data. This does not necessarily mean that these modes 
are spurious. As a consequence of splitting the data, the noise level 
increases, masking the low-amplitude modes. 

An amplitude change could be another reason for finding a 
mode in only one half of the data. In order to check this, we 
compared the semi-amplitudes of the modes derived independently 
from the two parts of the data, but with the frequencies fixed at the 
values given in TableQ The formal rms errors of semi-amplitudes 
derived from the time-series data similar to our dataset are known 
to be underestimated by a factor of about two (Handler et al. 2000, 
Jerzykiewicz et al. 2005; see also Montgomery & O'Donoghue 
1999). If we take this into account, denoting the doubled formal rms 
errors of the semi-amplitudes in a giver filter i by Si, we may con- 
clude that the differences in semi-amplitudes between the first and 
the second part of our data exceed 3 £j only for /1 in U (3.2 £u), V 
(4.7 e v ), and R (3.2 e a ), and for f 3 in U (4.0 eu) and V (3.5 ev). 
For about 70% of filter/mode combinations the differences were 
smaller than e,. We may therefore conclude that there are no sig- 
nificant amplitude differences between the two parts of the data 
except for a marginal change of the amplitude of two modes, /1 
(increase) and fz (decrease). 



4 THE AMPLITUDES 

A vast majority of the photometric observations of EC 14026 
stars was made without any filter. This was a compromise be- 
tween the required photometric accuracy and the exposure time 
needed for the very short periodicities observed in these variables. 
However, the multicolour observations are highly desirable since 
they can support mode identification, a starting point for success- 
ful asteroseismology. The best examples of the multicolour ob- 
servations of_EC 14026 stars made so far are the UBVR data 
of iKoenl ll998h for V2203 Cyg = KPD 2109+4401, the u'g'r' 
photometry of the same star and V 42 9 And = HS 0039+4302 = 
Balloon 8 4041013 Ijefferv et aljEooll) . four-band photometry of 
V 338 Ser iFalter et all2003l). and UB VR observations of LM Dra 
= PG 1618+563B ISilvotti et alfcOOch . 

Since Bal09 is relatively bright and mode identification is one 
of the main goals of our study, we observed it through wide-band 
filters. There is, however, a price. The residual standard deviation of 
our U BVR observations, a measure of typical accuracy of a single 
data point, is quite large, and amounts to 0.015, 0.011, 0.010, and 
0.010 mag, respectively (Table0. Fortunately, owing to the large 
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Figure 12. Relative amplitudes, Ai/Axs, i = {B, V, R}, for the five 
largest-amplitude modes detected in Bal09: fx, large filled circles, f%, f%, 
fn triplet, open circles, and /g, dots. Eiror bars were calculated from formal 
errors of semi-amplitudes given (Tabled. 

number of the observations we gathered, the detection threshold in 
Fourier periodograms is reasonably low, 0.8-1.4 mmag. 

The semi-amplitudes in all filters are given in Table □ It 
is worth no ting that, along with the main mode in V338 Ser 
(Koen et al. 1998a), the fx mode has the largest amplitude of all 
modes known in EC 14026 stars. Generally, the semi-amplitudes 
become smaller towards longer wavelengths as the theory predicts 
(see, e.g., Ramachandran et al. 2004). For five modes with the 
largest amplitudes we plot the amplitude ratios in Fig. 1121 They 
are normalized to the U -filter amplitude. 

The phase differences (Table are equal to zero within the 
errors for all modes. 



5 DISCUSSION 

5.1 Preliminary mode identification 

There are at least two possibilities to identify modes solely from the 
photometry. The first method relies on the comparison of the am- 
plitude ratios and phase differences with those calculated from the 
model atmospheres for a given mode. The amplitude ratios were al- 
ready used to constrain degree £ in EC 14026 stars by Koen 1 1998) 
and Jefferv et al. 1 2004). How ever, t hese papers and the theoretical 
work of Ramachandran et al. 1 2004) indicate that when going from 
shorter to longer wavelength in the optical domain, the amplitudes 
of low-degree (£ < 2) p modes for EC 14026 models behave simi- 
larly. Therefore, the £ for the I < 2 modes is hardly derivable from 
amplitude ratios. Modes with 1 = 3 and 4 are much better separated: 
these with £ = 3 show rather flat dependence on wavelength, while 
those with £ = 4, much steeper than the £ < 2 modes. The amplitude 
ratios for five modes shown in Fig. ll2l exhibit similar change with 
wavelength as the £ < 2 modes. We may therefore rather safely 
interpret them as modes having £ equal to 0, 1 or 2. 

The other possibility for mode identification from the photom- 
etry is the best match of the observed freq uencies to those calcu- 
lated from the model iBrassard et alJbOOll) is a good example]. We 
make here only a rough comparison of the observed periods with 
those calculated in the currently available models, giving some in- 
dications as to the possible identifications. A detailed mode identi- 
fication including the stability analysis of models in a wide range 
of input parameters needs much better modelling and is beyond the 
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Table 3. Parameters of interpolated models that have the period of the fun- 
damental radial mode equal to 356.19 s. 



Model 


C02 model 
sequence 


Age* 
[Myr] 


Evol. 
phase 


logT eff 


logs 


A 


6 


4.51 


HeCB 


4.403 


5.51 


B 


5 


24.81 


HeCB 


4.415 


5.51 


C 


4 


64.66 


HeCB 


4.436 


5.50 


D 


4 


107.02 


end of HeCB 


4.467 


5.51 


E 


5 


109.51 


end of HeCB 


4.475 


5.51 


F 


6 


104.51 


end of HeCB 


4.486 


5.51 


G 


7 


106.70 


CHeE 


4.499 


5.52 


* since ZAEHB. 



scope of this paper; we postpone such a thorough analysis to a sep- 
arate study. 

Let us start with the assumption that the main mode, fx, is the 
radial fundamental one. In addition to the expected wavelength- 
amplitude ration dependence, the justification for this assumption 
is its large amplitude. With this assumption, we can go to mod- 
els in order to try to interpret the other frequencies. A large set 
of pulsational properties of evolutionary models was published by 
Charpi net et alJ t2002) (hereafter C02) for seven extended horizon- 
tal branch (EHB) sequences. All sequences were calculated for a 
star with mass ~0.48 A/0, but with different envelope masses: 
0.0001, 0.0002, 0.0007, 0.0012, 0.0022, 0.0032, and 0.0042 M© 
for sequences 1 to 7, respectively. These models are shown in the 
log T e g - log g diagram in Fig.|13l Th e posi t ion o f Bal09, accord- 
ing to the parameters derived by Orei ro et all 120041) . is also shown. 

C02 provide also the pulsation periods for modes with £ up to 
3. The periods of the fundamental radial mode (£ = 0, po) are shown 
for all but the first two sequences in the bottom panel of Fig. fT3l The 
line for Pi = l//i = 356.19 s intersects the model values in seven 
points corresponding to seven models belonging to sequences 4 to 
7. Their properties, derived from linear interpolation of the model 
grid of C02, are listed in Table [3] We see that none of the models 
with envelope mass less than ~0.001 M© can be fit with the period 
of the fundamental radial mode equal to Pi. The seven models in 
the sequence from the least evolved to the most evolved we denote 
with the letters A to G. The A to C models cover the main part of the 
EHB helium core burning evolution (HeCB in Table|3j, D to F, its 
final stage, and G, the core helium exhaution (CHeE). The models 
are plotted in Fig. ll3l with filled gray circles. Note that regardless of 
the model, the theoret ical value of log g ~ 5.51 agrees reasonably 
well with the result of Oreiro et al. (2004). 

Having interpolated the models that match exactly Pi , we can 
try to compare the other periods as well. This is done in Fig. 1141 
for all seven models. For clarity, the £ = 3 modes are not shown in 
this figure. We see that the observed pattern resembles to some ex- 
tent the theoretical one (certainly with many high-overtone modes 
missing). We can therefore derive some conservative conclusions 
based on this comparison. Since besides the strong 356-s term, we 
have only three narrow ranges in period where the observed modes 
occur, we can try to constrain the possible identification. 

The separation between fx and the central term of the rotation- 
ally split triplet (f$) is best reproduced by model D. This is also the 
model with no £ < 2 modes having periods in the range between 
270 and 340 s, in accordance with observations. Moreover, model 
D is closest in T g to the observed value (Fig. 1131 . However, we 
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Figure 13. Top: Seven evolutionary sequences (labeled from 1 to 7) show- 
ing the extended horizontal branch evolution of a 0.47-0.48 Mq star, taken 
from|Charpinet et al. 1 2002). It is mainly the hydrogen envelope mass that 
differentiate the models; it ranges from 0.0001 Mq for sequence 1 up to 
0.0042 A/0 for sequence 7. The large dot with err or bars denotes the param- 
eters of Bal09 as derived by Oreiro et al. 1 2004). Seven filled gray circles 
show the position of models A to G for which the period of the fundamental 
radial mode is equal to Bottom: Periods of the fundamental radial 

mode, I = 0, po for five of the seven evolutionary sequences shown above. 
The dashed horizontal line corresponts to Pi = = 356.19 s. Gray filled 
circles show the positions of the seven models, A to G (see Table 13*1. that 
match the period exactly. 



are far from indicating that model D is the best one, and constrain- 
ing the evolutionary status of Bal09. The problem requires better 
modelling. 

Nevertheless, some indications as to the mode identification 
can be given even from this simple exercise. The mode closest in 
period to po, on the short-period side, is in all models the I = 1 
mode. It has radial order n = 1 for all models but G, where it has 
n = 2. However, while for the three less evolved models it is the pi 
mode, it is replaced in period for the remaining, more evolved mod- 
els, by (?i (D-F) or even 172 (G) due to the phenomenon of avoided 
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Figure 14. Comparison of the pulsation periods for models A-G, indicated 
in Fig. l 131 with periods detected in Bal09. For the models, the radial (I = 0) 
modes are shown with continuous lines, I = 1 modes with dashed lines, and 
I = 2 modes with dotted lines. See text for more details. 



Table 4. Identification of the family (p, f or g) and radial order (subscript) 
of modes with i = 0, 1 , and 2 closest to each of the three period groups 
detected in Bal09 for seven models listed in Tablel3l 



Period group 

350-356 s 



Period group 

263-265 s 



Period group 
214-215 s 
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pi 
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crossing. The g modes under consideration have mixed character 
and are typical for evolved models also in main-sequence pulsators. 
An £ = 1 identification can be therefore attributed with high prob- 
ability to the rotationally split modes /a, fa, and f^; the azimuthal 
orders would then be m = — 1, 0, and +1, respectively. 

There is a general trend that for the three groups of periods de- 
tected in Bal09, that is, around 350, 260 and 215 s, there are always 
nearby theoretical modes with £ = 0, 1 and 2. This can be seen in 
Fig-EI In Table|4] we identified the theoretical modes closest to the 
observed ones for all seven models. 



5.2 Coexistence of low and high-frequency modes 

There is no doubt that the low-frequency modes, presumably high- 
order g modes, exist in Bal09. The question arises whether such 
modes are also present in the other cool and evolved EC 14026 
stars. The low-frequency region is sometimes ignored in the fre- 
quency analysis and signals, even when found, are frequently at- 
tributed to the atmospheric transparency variations or removed dur- 
ing the detrending process. However, high-order g modes discov- 
ered so far have periods shorter than 2 hours. Thus, two to four 
cycles may be covered during a single night and these signals 
can be easily distinguished from the transparency variations, es- 
pecially in the differential photometry. It is therefore worthwhile 
to reanalyze some already existing long data sets with the aim 
of searching for low-frequency modes. While writing this paper, 
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we have learned about the discovery of a long-period mode in 
HS 0702+6043 ischuh etal1 l2004 . It is worth noting that both 
Bal09 and HS 0702+6043 are neighbours of the PG 1716+426 
group of pulsating subdwarfs in the log T c g - log g diagram, and 
some other EC 14026 stars are close by (V 338 Ser, V 391 Peg and 
KL UMa). These stars should be considered as primary targets in 
the search for high-order g modes in EC 14026 stars. 



5.3 Frequency splitting 

Modes close in frequency were detected in many EC 14026 stars, 
and some of them were supposed to be rotationally split ones. This 
would allow estimating rotational periods. However, up to now, 
the only rotationally split triplet was found in EO Cet = PB 8783 
(O'Donoghue et al. 1998). In several other EC 14026 stars only 
close doublets were found. The equidistant triplet found by us in 
Bal09 is therefore the second one detected in an EC 14026 star. If 
we accept that it is a rotationally split £ = 1 mode (see Sect. 5.1), we 
can translate the observed mean splitting, A/ = 1.58 /iHz, into the 
rotation period P rot = (1 — C„f)/A/, where C„e is the Ledoux 
constant. For all models considered in previous section but model 
D, the values of C n e provided by C02 lie between 0.01 and 0.04, 
implying P rot = 7.1 ± 0.1 d. Only for model D, the C'„e is higher 
and can reach 0.2. In that case P rot would be smaller, about 6 d. 
Bal09 is therefore very likely a slow rotator. The 7.1-day rotation 
period for a star with a mass of 0.48 Mq and log g = 5.51 results 
in the equatorial rotation velocity of 1.44 km/s. 



5.4 Final remarks 

It is clear for us that Bal09 is one of the most interesting targets for 
the future study by means of asteroseismology. The main reason 
for this is the coexistence of the pulsations inherent to both classes 
of pulsating sdB stars: the EC 14026 and PG 1716+426 stars. As 
far as the excited pulsations are concerned, it seems simpler than 
the best studied EC 14026-type star V 338 Ser that is a fast rotator 
( Hebe r et alll999h . Bal09 rotates slowly and is relatively bright. It 
is also exceptional because of the presence of a rotationally split 
mode and the fact that a large number of combination frequencies 
are detected. 

In addition to a better modelling we are going to undertake 
in near future, the next step in our understanding of Bal09 could 
be taken if a new photometric and spectroscopic campaign were 
organized. Such a campaign should at least remove the 1 d _1 am- 
biguity of some of the derived frequencies and lower the detec- 
tion threshold down to 0.1-0.2 mmag. Both these goals could be 
achieved if the star were observed in several observatories spread 
over the longitude, preferably in white light (or very wide band fil- 
ters). Multicolour observations can support mode identification, but 
we showed in Sect. 4 that the £ value can be reliably constrained in 
that way only for the largest-amplitude modes. It seems, anyway, 
that even with accurate amplitudes it will be hard to disentangle 
modes with t = 0, 1, and 2, and it might be necessary to go to 
satellite UV to make a significant progress (see Fontaine & Chayer 
2004). 

The campaign we plan should be supplemented by spectro- 
scopic observations. Some time-resolved spectroscopic observa- 
tions of EC 14026 stars have been already done. They were, how- 
ever, usually too short in duration or badly distributed in time to re- 
sult in a breakthrough in the mode identification. In case of Bal09, 
however, even confirming that the main frequency, /i , corresponds 



to the radial mode, as we suggested, would be helpful. The best 
way towards proper identification would come, however, from ac- 
curate matching the model frequencies to the observed ones. This is 
even possible with our results provided that a thorough modelling 
would be done. As a by-product, stellar parameters of Bal09 could 
be obtained. As noted above, we postpone this work to a separate 
paper. 
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